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FREE OSCILLATIONS OF AN ATMOSPHERE IN WHICH TEMPERATURE 
INCREASES LINEARLY WITH HEIGHT 
By C. L. Pekeris 


SUMMARY 


It is shown that when the temperature in the atmosphere increases 
linearly with height, the speed of propagation of long waves does not 
approach a limit with increasing wave length, as in the case of an 
atmosphere in which the temperature at great heights is assumed to he 
constant or decreasing, but increases linearly with the period. The 
group velocity ultimately also increases linearly with the period and 
becomes equal to half the phase velocity. The region of maxi mum energy 
of the oscillation is shifted to increasingly higher elevations as the 
period is increased. Whereas, in an atmosphere where the temperature 
gradient at great heights is negative or zero, the tides are similar to 
those in a uniform ocean of equivalent depth H (about 8 km, depending 
on the assumed vertical temperature distribution), the superposition of 
an outer envelope with a positive temperature gradient introduces a 
radical change into the nature of the tide. Insofar as it is still 
legitimate to refer the tides in such an atmosphere to those in an ocean 
of equivalent depth H, it can be said that H becomes a function of 
the period, which increases indefinitely with the period. The bearing 
of these results on the resonance theory of atmospheric tides is discussed. 


INTRODUCTION 


The fact that the solar semidiurnal tide in the atmosphere is about 
100 times larger than the equilibrium value led Kelvin to the conclusion 
that the atmosphere possesses a free period of tidal, oscillation of 
12 solar hours. The resonance period of the atmosphere must be within 
4 minutes of half a Bolar d&y in order to account for the observed ampli- 
fication and for the fact that the moon, whose tidal force is more than 
twice as large as that of the sun, excites a barometric oscillation having 
an amplitude of only one- sixteenth that of the solar wave. If the vertical 
temperature distribution in the atmosphere, assumed horizontally strati- 
fied, is known, it is possible to compute its free period of tidal oscil- 
lation following a method due to the work of G. I. Taylor (reference l). 
According to this method one first determines the speed of propagation of 
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long waves . V in a flat atmosphere having the same vertical temperature 
distribution. " The resonance period or periods of tidal oscillations of 
the atmosphere are then identical with those of an ocean of depth H 
enveloping the earth and are determined from the relation 

V = (gH) 1 / 2 (1) 

The dependence of the free period of the ocean on its depth H is 
known from classical tidal theory. For the solar semidiurnal oscilla- 
tion H = 7.87 kilometers ± 3 percent, and therefore V must he 
0.278 kilometer per second or 0.817 cqj where cq denotes the velocity 
of sound in the air next to the ground, that is, 3**0 meters per second 
for an assumed surface temperature of 288° K. On assuming a factor 
exp i(ot - kx) the divergence X of the wave motion is in the limit of 
long waves found to satisfy the equation 



c 2 = 7RT 
H = (^/gk 2 
V = cr/k 


( 2 ) 


subject to the boundary conditions that at the ground the vertical com- 
ponent of the velocity w vanishes; 


w = (c 2 Vg) + h[^ 7X + (c 2 /g)(dx/dzj] = 0 (3) 

and that the energy of wave motion per unit column of the atmosphere is 
not infinite. If p(z) denotes the density of the air as a function 
of height z then of the two independent solutions of equation (2) the 
one is accepted which gives a finite value for the energy integral U 

U = (l/2) J p(z)^u 2 + v 2 + w 2 ) dz (4) 


where u, v, and w denote the components of velocity due to the 
wave motion. 
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When the temperature at great heights is assumed to he sufficiently 
low, or decreasing, it is actually found (reference 2) that one of the 
solutions of equation (2) gives a finite value for U, while in the 
other the distribution of (u 2 + v 2 + w 2 ) with height is such that U 
becomes infinite as the integral is extended to infinity. When, however, 
the temperature at great heights is assumed to increase linearly without 
limit, it is found that for both solutions of equation (2) U becomes 
infinite, so that the question of the choice between the two solutions 
of equation ( 2 ) is left open. This is shown in a following section. 

• 

This work was conducted under the sponsorship and with the financial 
assistance of the National Advisory Committee for Aeronautics. 


SYMBOLS 


A constant in equation (63) 

A, J , B, Y coefficients in equations (8) and (9) 

b coefficient in equation (52) 

C = sin nb/sin jt(6 + n) 


c 0 
D = 

E = 
F 

F l> 

F 3> 

F(e) 

g 

H 

I 


velocity of sound at any level 
velocity of sound in air next to ground 
(2/7) x [(1 + njo] 1 / 2 

(8*/g)(l + ^oQ 1 + sOcd]" 1 / 2 

confluent hypergeometric function 
F 2 , G defined by equations ( 33 ) 

Fip N defined by equations (3^) 

Airy function (reference 3 ) 
acceleration of gravity 
depth in a uniform ocean 
Bessel function with Imaginary argument 
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J 


K = 

k 

M 


m = 


n = 


Bessel function with real argument 
F(l + 6) sin Jt(5 + n)/2 n+ ^- sin Jtn 

wave number 

coefficient in equation (47) 

[}.+ (7 - 1)(1 + n)]y7H 

(g/Rp) - 1 


P 

P B 

Q(u) 

R 

T 

t 

U 


u, v, w 
V 

W(x) 

x, z 
x 


*0 

y = 2 JEix. = 2^mz 


period of wave 

Brunt's period (reference 4) 

function defined by equation (5*0 

gas constant 

temperature 

timej in appendix B, the plane t 
energy integral 

components of velocity due to wave motion 
speed of propagation of long waves j phase velocity 
Whittaker's function (reference 5) 
defined by equations ( 19 ) and (20) 

distances in x- and z-directions 
roots of equation (35) 
value of x at ground 


z 0 


level above which temperature increases linearly* 
with height 


P 

r 


constant in equation (5) 
function in equation ( 19 ) 
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y ratio of specific heats 

5 = -1 - (n/2) + Jt (1 + n)/2r] + (os/^t) 

e function in equation (46) 


0 = (l/6)8 -1 / 2 |4S - 
P 

PO 

a 

t = a 2 /gk = kH 


xp / 2 

density of air 
undisturbed density 
frequency of wave 


•r^, Tg roots t of equation (l6) 

cp-L s 2 /mz - (nit/2) - (3«A) 
qpg = 2^mz + (nrt/2) + (rt/4) 


X 


divergence of wave motion 


ffl » 7 + (7 - 1)(1 + n) 


APPARENT FAILURE OF ENERGY- INTEGRAL CRITERION FOR AN ATMOSPHERE 
IN WHICH TEMPERATURE INCREASES LINEARLY WITH HEIGHT 


Let the temperature T increase linearly with height above the 
level- Zq, 


T = 0z, z > z 0 

c 2 = 7RT - 7gz/(l + n) 



n = (g/Rp) - 1 


(5) 


then one easily finds that the density decreases with height according 
to the power law 


p/p Q = (z/z 0 )“ n ~ 2 


(6) 
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Equation (2) then takes on the form 

(d 2 x/dz 2 ) - (n/z) (dx/dz) + (m/z)x = 0 (7) 

m = [7 + (r - 1)(1 + 

whose solution is 

X = Az( 1+n )/ 2 J 1+:a (2^) + Bz( 1+u )/ 2 Y 1+n (2\Jmz) (8) 



when n is an integer, or 

X = Az( 1+i1 V2j i+ii (2^) + (9) 


when n is not an integer. 


In evaluating the energy integral U in equation (4) for the- purpose 
of making the choice between the solutions in equation (8) or in equa- 
tion ( 9 ), one has, for a wave progressing in the x-direction. 


u = (i/gffl)(c 2 X - gv^! 
v = 0 


( 10 ) 


Since c 2 X ~ zX while, for large values of z, dx/dz s X it follows 
that at great heights both u and w vary like zX. Now for large 
values of z 


J^ +n (2\|mz) £ ^« 2 mz) - ^-A cos <p^ 
q> s [2\^z - (njt/2) - (3*A)J 

J .^(2^) ~ (n 2 mz)“^/ 1< ‘ COS <Pg 

cp 2 = J^2\jmz + (mt/2) + (*A)]] 


(H) 


( 12 ) 
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It follows that (u 2 + v 2 + w 2 ) in the integrand of equation (4) varies 
like z-*-/ 2 cos^cp^ or z^-/ 2 cos^cpg for the two solutions of equa- 
tion (9), respectively. The kinetic energy per unit volume therefore 
increases ultimately with height like z l/2 

and the energy integral U 

does not converge for either solution. . Hot only does the energy-integral 
criterion not enable one to make the choice between the two solutions of 
equation (2) but also neither solution yields a finite value for U. 


RESOLUTION. OF DIFFICULTY WITH THE ENERGY INTEGRAL 


The result obtained in the previous section that for an atmosphere 
in which the temperature increases linearly with height the energy inte- 
gral of the wave motion diverges for both solutions of the wave equa- 
tion is hard to reconcile with the following physical considerations. 

A medium such as the atmosphere, which is bounded by a perfectly- 
reflecting surface. at the ground but is exposed to space on the other 
side, would be expected to be facilitated in sustaining a free oscilla- 
tion the more it is capable of diverting wave energy from great heights 
toward the ground. The more the upper atmosphere behaves like a reflector, 
the greater will be its trapping power for waves, and, with it, its 
ability to propagate waves to great ranges - a characteristic of free 
oscillations. Now it is known (references 6 to 8) that an atmosphere of 
constant temperature throughout can completely trap waves either in the 
whole spectrum or in a limited frequency range. This is also true for 
an atmosphere of the Taylor type in which the temperature is assumed to 
decrease at a constant rate in the troposphere and to remain constant in 
the stratosphere. Now when the temperature is constant, the rays ema- 
nating from a source situated in the atmosphere are rectilinear, and all 
but the horizontal rays are lost to space. In the Taylor type atmosphere 
the rays are even bent upward in the troposphere, event ually to become 
rectilinear in the stratosphere. In spite of these nonconservative prop- 
erties, both types of atmosphere are oapable of sustaining free oscil- 
lations in certain frequency ranges. In the case of an atmosphere with 
a linearly increasing temperature all the rays are bent downward, so 
that its trapping powers should be even greater. Yet, the divergence of 
the energy integral found in the previous section indicates a tendency 
to allow the waves to leak: out to higher elevations. 


The resolution of this paradox lies in the fact that equation (2) 
for the divergence X is not exact but is an approximation for the 
limiting case of long waves. This approximation breaks down when the 
temperature increases indefinitely linearly with height. The exact 
equation for X is equation (13) 




r , x 

— 

(X as 2 . Zf) 

^+-i 

WfiS - 1] + 

1 1 dc 2 . g (7 - 1)1 

\c 2 az c 2 / 

dz 

' 

\c 2 / 

H c 2 dz + Hc 2 J 


•X= o -O'S) 
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(1 - H 2 k 2 )v = (c 2 x/g) + |h -7X + (c 2 / g) (dx/dz^J ^ (l 1 *-) 

replacing equation (3) • Now in the l im it of long wave s , when the wave 
number k vanishes, the long-wave approximation is made hy dropping 
the k 2 -term in equation (13). This approximation is valid in an atmos- 
phere of constant or decreasing temperature. In the case of an atmos- 
phere in which the temperature (and with it c 2 ) increases linearly with 
height, the terms in brackets in equation (13) decrease like l/z and 
eventually become smaller than the k 2 — term. In fact at extremely large 
heights equation (13) reduces, for a given k, to 

(d 2 x/dz 2 ) - lAt = 0 (15) 

indicating that the appropriate solution for X ultimately varies like 
e _kz not accor ding to the power law that would follow from equa- 

tions (9), (11), and (12). When the complete equation (13) is used one 
finds that the solution having the factor e -kz yields a finite value 
for the energy integral U in equation (4), whereas the other solution, 
with the factor e^ 2 , must be discarded because it gives a divergent 
value for U. With the complete wave equation (13) there is therefore 
no ambiguity as to the choice between the two independent solutions, nor 
is there any difficulty with the energy integral. 

Using the notation of equations (5) and putting 


x = 2kz 

t = 0 2 /gk = kH 

5 = -1 - (n/2) + [t( 1 + n)/2 + (u>/27t) 

<a = 7 + (7 - 1)(1 + n) J 

equation (13) reduces to 

x^X/dx 2 ) - n(dx/dx) + js + 1 + (n/2) - (x/4)]x = 0 


(16) 


(17) 


I 
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Choosing the solution of equation (17) which is houn de d at infinity, 

X = ^ /a »« (2B+2Hl) (1+n) {1 - [8< S * 1 + ”)/i + • • •} 

2 ’ 2 ( 18 ) 

where W(x) is Whittaker's function (reference 5* P- 337)- 53s® expo- 

nential behavior of X for large values of x assures the convergence 
of integral (4) for U, whereas in the second solution of equation ( 17 ), 

which behaves like x -1-5 e x / 2 at Infinity, the integral diverges. 


LIMITING POEM OF S0LDTI0N FOR LONG WAVES 

With the solution of the exact equation (13) given by equation ( 18 ) 
it becomes of interest to investigate the form It assumes in -the limit 
of long wavelengths. In particular the relation of this form 

of equation ( 18 ) to solution (9) of the approximate equation (2) is to 
be investigated, with a view of throwing light on the question of the 
choice of the proper linear combination of the two independent solutions 
appearing there. For this purpose the following relations (reference 5, 
PP- 337* 338, and 3^6) are used: 



where F denotes the confluent hypergeometric function. Thus, for the 
case when n is not an integer there is obtained 


X(x) = [r(-l - n)/P(-5 - n - ljJ^+^e-VapC-S, 2 + n, x) + 
|£(1 + n)/T(-6)]e-x/2F(-5 - n - 1 - n, x) 


( 21 ) 
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X(x) = |[P(2 + 8 + n) sin *(8 + + sil1 xr L+n e" x / 2 F(-8, 2 + 

n, x) - Ji(l + n) X I*(l + 8) sin (it8)/«je _x /^F(-8 - 1 - n, -n, x) 

( 22 ) 

-where in the last; -transformation "the following relation was used 


r(x) X I*(l - x) = Jt/sin 

« 

Passing now to the limit of long waves, 

t = kH « l' 

8 — »<o/27t = m/2k » 1 


(23) 


(24) 


With 8 very large, the F-functions in equation (22) reduce to Bessel 
functions: 

F(-8, 2 + n, x)— >|l - 8x/(2 + n) + (8x) 2 /2(2 + n)(3 + n) + • • -j = 

T(2 + n) ( 8x) -( n+1 )/ 2 J n+1 (2\f6x ) (25) 

F(-8 - 1 - n, -n, x)— »r(-n)(8x)( n+1 )/ 2 J_ 11 _ 1 (2\|8l) (2 6 ) 


It is important to observe that these limiting forms of the F's are 
valid only, when both 8 » 1 and x « 8, for otherwise the series 
for F depends on the higher powers of x, when it iB no longer legiti- 
mate to replace the factor 8 - k by 8 in the coefficients. Sub- 
stituting now equations ( 25 ) and ( 26 ) into equation (22) and performing 
some further reductions with the aid of equation ( 23 ), there is obtained 


X = K exp J^(y^/8S)^j [y I 1 + 1 J n +i(y) + Cy n+1 j_ n _ 1 (y)j 


(27) 
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y = 2^3x = 2.yjm, 

K = r(l + 8) sin Jt(8 + n)/2 n+ ^ sin Jtn 
C = sin «8/sin n(8 + n) 


(23) 

(29) 


Except for the factor exp [j(y^/8&)] , which is slowly variable in the 
limit of large value? of 8, equation (27) is identical with equation (9)* 
The ratio of constants B/A in the latter is now> however, no longer 
arbitrary but is a definite function C of 8, as giveh in equation (29) . 

The goal of determining the appropriate solution of the wave equation 
for an atmosphere with a linear increase of temperature has thus been 
achieved, but on closer study the physical implications of the results 
are rather disturbing. In order to ascertain the physical significance 
of the parameter 8 the phase velocity and the period are determined. 

The phase velocity V is given by 


V/< c 0 = cr/kc 0 = (2/7) x [(I + n)cu] 1 ^ 2 (l/y) s D/y (30) 

while the period of the wave P, measured in' seconds, is given by 

P = 2j r/cr = (8*/g)(l + n)c 0 [Tl + n)tu] ^(b/y) e E8/y (31) 


For n = k.5, then D = 6.36 and E = 1080, while for n = 8 . 5 , which 
is appropriate for the positive lapse rate prevailing in the NACA model 
atmosphere between 83 and 120 kilometers (reference 9 ), D = 10.0 and 
E = Il80. For v/cq of the order of unity, the period P changes by 
about 3 minutes as 8 changes by unity. But, since C(8) in equation ( 29 ) 
is a periodic function of 8 of period unity and takes on all values 
between -00 and <» as 8 changes by 1, it follows that in the li m it of 
long periods X does not approach a fixed form but changes its character 
radically as the period is changed by only a few minutes. This increased 
sensitivity of X to the period for an atmosphere with positive lapse 
rate is to be contrasted with the decreasing sensitivity to changes in 
period for atmospheres with zero or negative lapse rates. 
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FRR'K OSCIIIATIONS OF AH ATMOSPHERE IH WHICH TEMPERATURE 
TW CRRARER T.TWRfl 'RT.T WITH HEIGHT 


In order to gain more insight into the acoustic propagating prop- 
erties of an a tmo sphere with a constant positive temperature gradient, 
the free oscillations or such an atmosphere when it extends from the 
ground up will he investigated. Therefore the variation of the phase 
velocity V with the period P for the first and far the higher normal 
modes will te studied. Such a variation is shown, far example, hy 
curve I in figure 2, for the first mode of an atmosphere with constant 
negative temperature gradient (references 8 and 10) . For a Taylor type 
atmosphere having the same negative temperature gradient in the tropo- 
sphere, the corresponding curve starts with a value of 0.86 at 
P ss 2 minutes and reaches its limiting value of about 0.92 at 
P = 4 minutes (references 8 and 10) . 

The divergence X from equation (22) can he written in the form 


X = - Ge _x /% , 2 


(32) 


with 

G a |r(l-t-ii)r(2+n)P(l + 6) sin JCnXsin «B]/[*r(2 + 5 + n) x Bin it(8 +nj] 

F^ = F(-8, 2 n, x) 

F 2 = F(-8 - 1 - n, -n, x) 

(33) 

The possible free oscillations of the atmosphere are determined from the 
condition, of the vanishing of the vertical component of velocity w at 
the ground., which, according to equation (l4), is given hy 


(1 * m)(l - t 2 )(w/ 7HX) = -1 - n + (x/2r) + (x/x)(dx/dx) = 0 (34) 


With 

F 3 = F(1 - 5, 3 + n, x) 

Fij. = F(-5 - n, 1 - n, x) 

D = - Gjr 1-n F 2 

N = [(1 + nj/xj]?! - [i/(2 + n[[F3 - G [(8 + 1 + nj/njx- 1 " 1 ^!,. 


(34a) 
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equation (34) leads to the condition 


(1 + n)(l - T^w/yHX) = -(1 + n) + (x/2) J(1/t) - l] + xH/D = 0 
x = *0 


( 35 ) 


X = Xq 

where xq denotes the value of x at the ground. 


The problem of determining the dependence of the phase velocity V 
on the period P involves the determination of the dependence of the 
roots x of equation (35) on the parameter 5. For a given 6 one 
■first determines the two roots f of equation (l6): 


T i = ^7[8 + 1 + (n/2)] + { r2 E + 1 + (n/2)] 2 - (1 + n -Da + n) 

t 2 = + 1 + (n/2[| - Jr 2 [8 + 1 + (n/2)] 2 - (1 + n)o>j- ^^1 + n) 

(37) 

For each r one finds the roots x of equation (35) and determines V 
and P from 


1/2 


V/c 0 = [i(l + n)T /7x] 


P = (2nc 0 /g) [i(l + n)/7Tx] 


1/2 


(38) 

(39) 


The free oscillations of an atmosphere have been determined with 
n = 4.5, corresponding to a positive temperature gradient of about 
6. 4° 'per kilometer. Figure 1 shows the variation of the roots x with 
8 for the two branches and t 2 . It is shown in appendix A that 

for the first mode in the_ "J^-branch, which is shown by curve A in fig- 
ure 1, the lower end of the curve is given by 

x = 1.88(1 - 8) 1 / 6 * 5 (40) 

for sufficiently small values of x. With tq — 4/ll, it follows from 
equations (38) and (39) that along this end of curve A, P— >», 
v/cq— and V/cq— >(l/lO)P,_ when P is measured in minutes. Hence 
in an atmosphere with constant positive temperature gradient, the phase 
velocity does not approach a limit with increasing period but increases 
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linear ly with the period. The same result applies to the higher modes. 
Figure 2 shows the function V(P) which results from the x(5)-curves in 
figure 1. It is shown in appendix B that in the Tg-branch no waves can 
he propagated for periods less than 4.5 minutes, which is Brunt's period 
(reference *0 determined from 

i 

c 2 o 2 = g 2 (7 - 1) + g(dc 2 /dz) (4l) 

In all the modes v/c q grows indefinitely with the period. 

It is of interest to see the type of dispersion one gets by using 
the "long wave" solution (27) • Dropping the exponential term. 


X = y n+1 J n+1 (y). + Cy n+1 J_ n _ 1 (y) (1+2) 


dx/dy = y n+1 J n (y) - Cy^J.^y) (1+3) 

while the secular equation (35) reduces to 


-(1 + n) + (n^/lto) + (y/2) Jj n (y) - CJ_ n (yf|/[j n+ i(y) + CJ_ I1 _i(y)J = 0 

(41+) 

In this case one can solve directly for C, since all the other terms 
do not depend explicitly on 6: 


- " [T 1 + n)(2/y) r ( 7 y/ 2 n^ J n+1 (y)^J^|j_ n (y) + [(l + n)(2/y) - 

(7y/acD)Jj_ n _ 1 (y)J- (1+5) 


With C(6) given in equation (29), equation (45) allows one to obtain B 
as a function of y, from which one determines the function V(P) by 
losing equations (30) and (31). The results are shown in figure 3* 


It is seen that the general shape of the V(P) -curves is similar to 
that of the exact curves given in figure 2. The deviations are due to 
the fact that, in the upper right-hand corner, S is of the order of 
unity, for which the long-wave approximation cannot apply. The curves 
in figure 3 become also inaccurate as the cut-off period is approached. 
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because there y is of the order of &, anfl the Bessel function approxi- 
mation in relations ( 25 ) and (2 6 ) breaks down. The analysis for this 
limiting case is given in appendix B. 


APPLICATION TO THEORY OF ATMOSPHERIC TIDES 


This investigation has shown that in an atmosphere where the temper- 
ature increases indefinitely with height at a constant rate, the acoustic 
propagation properties are radically different from those that have 
hitherto been encountered in the theory of atmospheric tides. The 
principal new feature is that the phase velocity V no longer reaches 
a limit for waves of long period. For atmospheres where V reaches a 
limit, the theory of atmospheric tides is simplified by an application 
of Taylor's theorem (reference 1), according to which the tides in the 
atmosphere are equivalent to those in an ocean of uniform depth H = V^/g. 
In our case, V grows indefinitely with the period, and, insofar as it 
is still legitimate to refer to an equivalent ocean of depth H, one 
may say that H depends on the period, ultimately increasing as the 
square of the period. What happens is that, as the period is increased, 
the distribution of energy of wave motion with height no longer approaches 
a fixed form, but the energy continues to spread to higher elevations as 
the period is increased. A discussion will not be given here of the 
nature of the tides in an atmosphere with constant positive temperature 
gradient extending from the ground up. In such an atmosphere the tides 
are probably different from those in an ocean of fixed depth. An investi- 
gation of this problem would be in order if the evidence would continue 
to accumulate to the effect that the temperature in the ionosphere 
increases continuously with height. 

Pending such an investigation, suffice it to remark that there is 
no difficulty in treating the case where a layer with positive tempera- 
ture gradient is topped by an outer envelope of constant or decreasing 
temperature. In that case the proper solution for the layer is given 
by equation ( 9 ), with B/A arbitrary. It would also follow that in the 
absence of proof as to the temperature of the outer envelope of the atmos- 
phere one should not in theoretical studies allow an- assumed positive 
gradient in an interior layer to extend through the outer envelope. The 
consequences of such an extrapolation of the temperature curve, if rigor- 
ously determined, are, in the light of the results of this investigation, 

’ likely to affect seriously the nature of the atmospheric tides. Any 
preliminary results obtained for an atmosphere with a positive tempera- 
ture gradient extending through the envelope should therefore be checked 
by comparing them with those for a modified atmosphere in which, say, 
the temperature in the outer envelope is assumed to be constant or 
decreasing with height. It would also be helpful to compute in each 
case the distribution with height of the energy density of the wave 
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motion (p/2)(u 2 + v 2 + v 2 ) in order to assure that any segment of the 
temperature curve which is as yet not supported by observations does not 
effect a marked redistribution in the wave energy. 


Institute for Advanced Study 

Princeton, N. J., September 26, 19^9 
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APPENDIX A 

DEPENDENCE OF PEASE VELOCITY V ON PERIOD' P WHEN 
PERIOD IS INCREASED INDEFINITELY 


Only the t 2 -branch of the first mode will be treated here and in the ( 

particular case when n = 4.5*' The higher-order inodes, as well as those 
of the T}_-branch, can be dealt with similarly. The problem is to determine 
the limiting form of the root x(5) of equation (35) as 5—^1 - e. Then 


T 2 = k/l± 1 
G _» _ er(5.5)/6.5 

J 


(46) 


Also for a vanishing x all the F*s reduce to unity. Now assume that 


G— »Mx 6 -5 (47) 

as x— >0, and seek to determine the coefficient M from equation (35) . 
On retaining only leading terms, equation (35) becomes 


-5-5 + 0.875x + [p.5 - (x/6.5)](l + Mx) = 0 
On annulling the coefficient of x there is obtained 

M = -0.721 = -(1 - 8)[t(5.5)/6.5| x" 6 ' 5 

1/6.5 


m 


(49) 


x = 1.88(1 - 8) 


(50) 
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APPENDIX B 

DETERMINATION OF CUT-OFF PERIOD 


If equation (45) is evaluated for large values of y "by UBing the 
asymptotic forms of the Bessel functions it is found that all modes 
coalesce at a limiting period (5-7 min for n = 4.5), for vhich the 
phase velocity vanishes. This result can, however, not he relied upon 
since when y becomes of the order of 5 the Bessel function approxi- 
mations in relations .(25) and (26) are no longer valid. In order to 
determine the cut-off period the a s y mp totic form of W in equation (18) 
has to be obtained when both & and x become large and are of the . 
same order of magnitude. For this purpose Whittaker 1 s (reference 5> 
p. 339) countour- integral representation for W is used, which can 
be put in the form 

X s x n /2w(x) 

- (l/2xi)e ilt5 r(l + 5 ) e~ x / 2 x 1+n <^) f S" 1 ( 1 + t) 5+1+I1 e- ;)cfc at ( 51 ) 

the path starting at ~ and returning to « after encircling the origin 
in the positive direction, but excluding the point t = -1. Now put 


x = 45 - bS 
t = u - (1/2) 
then near t = -( 1 / 2 ) there is the expansion 


(52) 


r 8 -^! + t) 84,141 ^-* 6 = ( u + I) yfu - 0 exp [Ii*5 + (x/2) + Q(uj] ( 53 ) 


Q(u) = (Sb)u + ( 165 / 3 )u3 + (6 Us/5)u5 + . . . 


X = (l/2xi)r(l + Sjx 1 * 


n o| + 1) 1 - 1) 

U — 


exp Q(u) du 


(55) 


where the path in the u-plane is similar to the former path in the .t-plane. 
For large values of 5 and for values of b of the order of 5"^/3 , 
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tlie principal contribution to the integral in equation (55) comes from 
the immediate vicinity of the origin. This can be made evident by putting 


(l6s/3)u3 = v3 
A = (3/l6) 1 / 3 b& 2 / 3 


(56) 


when the u^-term in equation (5^) becomes of the order of 6“^/3, while 
the expression in brackets in the integrand of equation (55) reduces to 
-2“ n to within terms of the order of 6-1/3. 

Hence, 


X « -(l/2rti)2“ n x 1+n r(l + S)(3/l6B) 1 /3^ exp (v3 + Av) dv (57) 

Here the path starts at ooelV3 and terminates at °° after 

encircling the origin in the positive direction. The integral can be 
evaluated in terms of the Airy function (reference 3 , P- 189), yielding 



(58) 


(59) 


(60) 
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Equation (60) shows that heyond 4S, X(x) is no longer oscillatory. 
The largest zeros of X are determined from the zeros of F(0) (refer- 
ence 3, p. 751) in equation (59): 


0-L = 2.383 


0 2 = 5.610 


0 3 = 8.657 




(61) 


On substituting expression (59) into the secular equation (35) one finds 
that in the limit of large values of 5 the roots of equation (35) 
coincide with the roots of x- Hence the largest roots of equation (35) 
are given hy 

l/3 

x ± = 4s - (36S0i 2 ) ( g 2) 


where the 0-j/s are given in equation (6l) . From equation (62) it 
follows that near the cut-off period 


V/c 0 -^(A/6) [l + (l/8)(60i/B) 2 /3j 

P = p b [l + (l/8)(60i/B) 2 /3] 


(63) 


where A is a constant and Pg denotes Brunt 1 b period (reference 4) . 
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